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We study a system consisting of a junction of JV quantum wires, where the junction is characterized 
by a scalar S'-matrix, and an impurity spin is coupled to the electrons close to the junction. The 
wires are modeled as weakly interacting Tomonaga-Luttinger liquids. We derive the renormalization 
group equations for the Kondo couplings of the spin to the electronic modes on different wires, and 
analyze the renormalization group flows and fixed points for different values of the initial Kondo 
couplings and of the junction S'-matrix (such as the decoupled S- matrix and the Griffiths S-matrix). 
We generally find that the Kondo couplings flow towards large and antiferromagnetic values in one 
of two possible ways. For the Griffiths S-matrix, we study one of the strong coupling flows by a 
perturbative expansion in the inverse of the Kondo coupling; we find that at large distances, the 
system approaches the ferromagnetic fixed point of the decoupled S-matrix. For the decoupled S- 
matrix with antiferromagnetic Kondo couplings and weak inter-electron interactions, the flows are 
to one of two strong coupling fixed points in which all the channels are strongly coupled to each 
other through the impurity spin. But strong inter-electron interactions, with K p < N/(N + 2), 
stabilize a multi-channel fixed point in which the coupling between different channels goes to zero. 
We have also studied the temperature dependence of the conductance at the decoupled and Griffiths 
S-matrices. 

PACS numbers: 73.63.Nm, 72.15.Qm, 73.23.-b, 71.10.Pm 



I. INTRODUCTION 

The area of molecular electronics has grown tremen- 
dously in recent years as a result of the drive towards 
smaller and smaller electronic devices jjj . Molecular elec- 
tronic circuits typically need multi-probe junctions. The 
first experimental growths of three-terminal nanotube 
junctions were not well controlled 2]; more recently, new 
growth methods have been developed whereby uniform 
K-junctions have been produced 0,0,01 • Transport mea- 
surements have also been carried out for the Y- junctions 
0, as well as for three-terminal junctions obtained by 
mergingtogether single-walled nanotubes by molecular 
linkers pj. 

On the theoretical side, there have been several stud- 
ies of junctions of quantum wires. There have been de- 
tailed studies of carbon nanotubes with different pro- 
posed structures for the junction 0,0. Several groups 
have analyzed the geometry and stability of the junctions 
[lOlllll . Junctions of quantum wires have also been stud- 
ied E3, 0, EH IH, Ell m terms of one-dimensional 
wires, with the junction being modeled by a scattering 
matrix S. These studies include the effects of electron- 
electron interactions which are often cast in the language 
of Tomonaga-Luttinger liquid (TLL) theory [11 El Eg . 

Many earlier studies of junctions have only included 
'scalar' scatterings at the junction, i.e., the .S'-matrix 
has been taken to be spin- independent. The response 
of a junction of quantum wires to a magnetic impu- 
rity or an impurity spin at the junction has recently 
been studied both experimentally |2l| and theoretically 



HHH IH 0. As is well-known in three dimen- 
sions, an impurity spin can lead to the Kondo effect |27| . 
[The coupling between the conduction electrons and the 
impurity spin grows as one goes to lower temperatures; 
this leads to a larger scattering and therefore a larger 
resistance as long as the temperature is higher than the 
Kondo temperature Tr ■ Below Tr , the resistance due to 
scattering from the impurity spin decreases (if the value 
of the impurity spin S is larger than or equal to half the 
number of channels N) because the spin decouples from 
the electrons.] The Kondo effect for a 'two-wire junc- 
tion' in a TLL wire has been studied by several groups 

mmmmmilimm Using 

a renormalization 

group (RG) analysis for weak potential scattering, Fu- 
rusaki and Nagaosa showed that for an impurity spin of 
1/2, there is a stable strong coupling fixed point (FP) 
consisting of two semi-infinite uncoupled TLL wires and 
a spin singlet |2^. For strong potential scattering, the 
above FP is reached when the inter-electron interactions 
are weak. However, sufficiently strong inter-electron in- 
teractions are known to stabilize the two-channel Kondo 
FP instead [30|. The Kondo effect has also been stud- 
ied in crossed TLL wires [3fJ and in multi-wire systems 

mm 

In this paper, we consider a junction of quantum wires 
which is characterized by an S'-matrix at the junction; 
further, an impurity spin is coupled to the electrons at the 
junction. The wires are modeled as semi-infinite TLLs. 
The details of the model defined in the continuum will be 
described in Sec. II. In Sec. Ill, we will discuss how RG 
equations for the Kondo couplings and for the S'-matrix 
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at the junction can be obtained by successively integrat- 
ing out the electronic modes far from the Fermi energy. 
We find that the flow of the Kondo couplings involve the 
S'-matrix elements, but the flow of the S-matrix elements 
do not involve the Kondo couplings (up to second order 
in the latter). To simplify our analysis, therefore, we con- 
centrate on the FPs of the S-matrix RG equations and 
study how the Kondo couplings evolve in Sec. IV. For the 
case of N decoupled wires, we find that for a large range 
of initial values of the Kondo couplings, the system flows 
to a multi-channel ferromagnetic (FM) FP lying at zero 
coupling. This FP is associated with spin-flip scatterings 
of the electrons from the impurity spin whose tempera- 
ture dependence will be discussed. Outside this range, 
the flow is towards a strong antiferromagnetic (AFM) 
coupling. On the other hand, at the Griffiths S-matrix 
(defined below), there is no stable FP for finite values 
of the Kondo couplings, and the system flows towards 
strong AFM coupling in two possible ways. We also con- 
sider the case when the scattering matrix has a chiral 
form. In this case, we find that the Kondo coupling ma- 
trix for the three wire case has three independent degrees 
of freedom and a single FP at strong coupling. 

The strong coupling flows will be further discussed in 
Sec. V where we will consider some lattice models at the 
microscopic length scale. As in the three-dimensional 
Kondo problem, we find that there are various possibili- 
ties depending on the number of wires N and the spin S 
of the impurity, such as the under-screened, over-screened 
and exactly screened cases [39] . We will generally see 
that a Kondo coupling which is small at high tempera- 
tures (small length scales) can become large at low tem- 
peratures (large length scales). In Sec. VI, we will show 
that the vicinity of one of the strong coupling FPs can be 
studied through an expansion in the inverse of the cou- 
pling; we will then find that the large coupling can be 
re-interpreted as a small coupling in a different model. 

In Sec. VII, we will study the case of decoupled wires 
with strong interactions using the technique of bosoniza- 
tion. Analogous to the results of [^j, we nn( i that the 
multi-channel (N > 2) AFM Kondo FP is stabilized for 
K p < N/(N + 2). We will discuss the temperature de- 
pendence of the conductance in Sec. VIII at both high 
and low temperature; we will compare the behaviors of 
Fermi liquids and TLLs. Sec. IX will contain some con- 
cluding remarks. A condensed version of some parts of 
this paper has appeared elsewhere |2fjj . 

Before proceeding further, we would like to emphasize 
that we have not used bosonization in this paper (except 
in Sec. VII), although this is a pow erful and commonly 
used method for studying TLLs [13, |2Qj . In the pres- 
ence of a junction with a general scattering matrix, it is 
not known whether the idea of bosonization can be imple- 
mented. (Some reasons for the difficulty in bosonizing are 
explained in Ref. [13 )• It is therefore necessary to work 
directly in the fermionic language. We have adopted the 



following point of view in this work |l4 . |40| . We start 
with non-interacting electrons for which the scattering 
matrix approach and the Landauer formalism for study- 
ing electronic transport 0, 0] are justified. We then 
assume that the interactions between the electrons are 
weak, and treat the interactions to first order in pertur- 
bation theory to derive the RG equations. This is the 
approach used in most of this paper. Only in Sec. VII 
do we use bosonization to discuss the effect of strong in- 
teractions for the case of decoupled wires, since that is 
one of the cases where bosonization can be used. 

II. MODEL FOR SEVERAL WIRES COUPLED 
TO AN IMPURITY SPIN 

We begin with TV semi-infinite quantum wires which 
meet at one site which is the junction; on each wire, the 
spatial coordinate x will be taken to increase from zero at 
the junction to oo as we move far away from the junction. 

The incoming and outgoing fields are related by an 
S'-matrix at the junction, which is an N x N unitary 
matrix whose explicit values depend on the details of the 
junction. Hence the wave function corresponding to an 
electron with spin a (a =|, !) and wave number k which 
is incoming in wire i (i = 1, 2, • • • , N) is given by 

tpiak(x) = e - l{k+kF)x + S u e l{k+kF)x on wire i, 

= Sji e t{k+kF)x on wire j^i . (1) 

Here k is the wave number with respect to the Fermi 
wave number Uf (i.e., k = implies that the energy 
of the electron is equal to the Fermi energy Ep). Wc 
will take k to go from —A to A, where A is a cut-off of 
the order of kp ; we will eventually only be interested in 
the long wavelength modes with |fc| <C A. We will use 
a linearized approximation for the dispersion relation so 
that the energy of an electron with wave number k is 
given by vpk with respect to the Fermi energy; here vp is 
the Fermi velocity, and we are setting h = 1 . In Eq. Q , 
we will refer to the waves going as e~ lkx as the incoming 
part tpnak, and the waves going as e lkx as the outgoing 
part ipoiak or ipojak- 

The second quantized annihilation operator corre- 
sponding to the wave function in Eq. is given by 

^iak(x) = C iak 1pi a k(x) , (2) 

where the wire index i runs from 1 to A, and the total 
second quantized operator is given by 

& a (x) = [ ^ c iak ip lak (x) . (3) 

(Note that it is not possible to quantize the system in 
terms of N independent fields on each of the wires, be- 
cause an electron that is incoming on one wire has out- 
going components on all the other wires as well). The 
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non-interacting part of the Hamiltonian for the electrons 
is then given by 

H = v F J2 j ^ k c Lfe c »fe ■ ( 4 ) 

If the impurity spin is coupled to the electrons at the 
junction, that part of the Hamiltonian is given by 

H spin = ]T J S-Vt(x = 0) ?f-* fi (x = 0) , (5) 

where a denotes the Pauli matrices. (For simplicity, we 
will assume an isotropic spin coupling J x = J y = J z ). Eq. 
can be written in terms of second quantized operators 



H. 



id a,/3 



dk\ dki 
~2tt ~2tT 



where J io = J(l + J2i s u)( l + Em 5 "']) is a Hermi- 
tian matrix. In general, however, the impurity spin may 
also be coupled to the electrons at other sites which are 
slightly away from the junction; for instance, this may be 
true if the model is defined on a lattice at the microscopic 
scale as we will see in Sec. V. It is therefore convenient 
to take Jij to be an arbitrary Hermitian matrix which is 
not necessarily related to the entries of the S'-matrix in 
any simple way. 

Next, let us consider density-density interactions be- 
tween the electrons in each wire of the form (we will drop 
the wire index i for the moment) 



dx dy p(x) U(x - y) p(y) , (7) 



where the density p is given in terms of the second quan- 
tized electron field ^ a (x) as p = fy^^ + ^j^j.. As men- 
tioned earlier for the wave- functions, the electron field 
can also be written in terms of outgoing and incoming 
fields as 



(8) 



If the range of the interaction U (x) is short (of the or- 
der of the Fermi wavelength 2-K/kp), such as that of a 
screened Coulomb repulsion, the Hamiltonian in Q can 
be written as 

Hint = 



On 



<x,0 



(9) 



where 



fl i = U(2k F ) , 
and 52 = 9i = U(Q) 



(10) 



For repulsive and attractive interactions, gi > and < 
respectively. (We have ignored umklapp scattering terms 
here; they only arise if the model is defined on a lattice 
and we are at half- filling). 



III. THE RENORMALIZATION GROUP 
EQUATIONS 

It is known that the interaction parameters g\ , g2 and 
<?4 satisfy some RG equations 0|; the solutions of the 
lowest order RG equations are given by 
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(11) 



where L denotes the length scale. 

In general, the couplings <?i, 172 and 9^ can have dif- 
ferent values on different wires; hence we have to add a 
subscript i to them. For weak interactions, i.e., when gu, 
<72i and gu are all much less than 2ttv f , we can derive the 
RG equations for the S-matrix at the junction [l4l |40|. 
Let us define a parameter 



a, 



.92/ 



2 gu 



2ttv f 



(12) 



which is a function of length scale due to Eqs. (|11|) . and 
a diagonal matrix M whose entries are given by 



1 



(13) 



Then the RG equations can be written in the matrix form 
dS 



d\nL 



M - SM^S . 



(14) 



The FPs of this equation are given by the condition M = 
SM^S. 

We use the technique of 'poor man's RG' [liJEl to de- 
rive the renormalization of the S'-matrix and the Kondo 
coupling matrix J^. Briefly, this involves using the sec- 
ond order perturbation expression for the low energy ef- 
fective Hamiltonian, 



H t , 



E 



\h >< l 2 \H pelt \h >< h\H pC rt\h >< h\ 



Ei -E h 



(15) 
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J 



J 




7T S(k± 



hi) 



(17) 



(a) 



FIG. 1: Pictures of the terms which contribute to the renor- 
malization of the Kondo coupling matrix J to order J 2 and 
<?2 J respectively; <?2 denotes the coefficient of the electron- 
electron interaction. Thin lines and thick lines denote low 
energy and high energy electrons respectively. 



where the perturbation H pcrt is given by the sum of H spm 
and H[ nt in Eqs. © an d ©, h and h denote two energy 
states, and h denotes high energy states. We now restrict 
the sum over h in Eq. (|15f) to run over states for which 
the energy difference Eh ~ E\ lies within an energy shell E 
and E+dE; we have assumed that the difference between 
different low energy states is much smaller than E, so 
that we can simply write Ei t — Ei in the denominator 
of the above equation. We then see that the change in 
the effective Hamiltonian dH c g is proportional to dE/E 
which is equal to — dlnL, where the length scale L is 
inversely related to the energy scale E. In this way, we 
get an RG equation for the derivatives with respect to 
InL of various parameters appearing in the low energy 
Hamiltonian. 

Using this method, we find that the Kondo couplings 
Jij do not contribute to the renormalization of the S- 
matrix in Eq. Ijl4|) up to second order in Jij. (This is 
not true beyond second order; however, we will only work 
to second order here assuming that the Jij are small). 
On the other hand, the S-matrix does contribute to the 
renormalization of the Jij through the interaction Hamil- 
tonian in Eq. ©; this is because the relation between 
the outgoing field on wire i (i.e., ^oia) and the opera- 
tors Cj a involves the S-matrix. For instance, the terms 
involving g2i in Eq. © take the form 

w f A f A f A f A dkl dk2 dks dki 

http J —A J-A J-A J-A 2^" ~2T ~2T 
x 7T 5{kt - k 2 + k 3 ~ k 4 ) 92i 

X S *j C ]ak x 4(3k 2 C ipk 3 SuC lak4 , (16) 

where we have used the identity 



I 

Jo 



k\ — k 2 + k$ — &4 — ie 
1 



with e being an infinitesimal positive number. [In Eq. 
(|16f) . we have kept only the <5-function term and have 
dropped the principal part term since the latter can be 
either positive or negative, and its contribution vanishes 
when one integrates over the variables fc.j .] Note that the 
terms involving g 2 in Eq. (|l(jfl (as well as those involving 
gi and 54 in Eq. ©) conserve momentum while the 
Kondo coupling terms in Eq. © do not. 

We will omit the details of the RG calculations here 
apart from making a few comments below. We find that 



dJi : 
dh^L 



— [ J2 J * kJk i 



i P 



ki — + ks - k 4 



2nvp 

+ g 92i Sij ^ JikS* k + 2 92j Sjt ^ JkjSjk 
k k 

- 2 X] (- 92fe _ 2gik) (JikS£ k S k] + SfoSkkJkj) ] , 

k 

(18) 

where Sij is the S'-matrix at the length scale L. Eq. Ijl8|l 
is the key result of this paper. Note that it maintains 
the hcrmiticity of the matrix Jij . Eq. I|18|) always has a 
trivial FP at = 0. 

Tet us briefly comment on the origin of the various 
terms on the right hand side of Eq. (|18|l . The first and 
second lines arise from Figs. 1 (a) and (b) respectively. 
(The terms of order J 2 in the first line have been de- 
rived in Ref. [13] )■ The parameters gu and g^ do not 
appear in the second line of Eq. (|18fl since the terms 
which are proportional to these parameters either do not 
appear in the numerator of Eq. (|15[) because they are 
not allowed by momentum conservation, or they appear 
in Eq. (|15fl but their contribution vanishes because the 
Pauli matrices are traceless. Finally, the third line of Eq. 
(|T%|| arises as follows. In Ref. 0, the RG equation for 
the S'-matrix was derived. This was based on the idea 
that due to reflections at the junction (these arise from 
the diagonal elements of the S-matrix which are the re- 
flection amplitudes), there are Friedel oscillations in the 
density of the electrons; the amplitudes of these oscilla- 
tions are proportional to Skk and S£ k in wire k. We now 
treat the interactions in the Hartree-Fock approximation 
|l4|: this results in reflections from the Friedel oscilla- 
tions with a strength proportional to g 2 k ~ 2<?ifc in wire 
k. Now, an electron going from wire j to i can either (i) 
first go from wire j to wire k with a transmission ampli- 
tude Skj, scatter from the Friedel oscillations in wire k 
with an amplitude {gik — ^9ik)S kk , and finally scatter off 
the impurity spin from wire k to wire i with amplitude 
Jik, or (ii) first scatter off the impurity spin from wire 
j to wire k with amplitude J^j , scatter from the Friedel 
oscillations in wire k with an amplitude (<?2fc — 2gik)Skk, 
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and finally scatter from wire k to wire i with a transmis- 
sion amplitude <SL. These two processes give rise to the 
third line of Eq. ftty 

It is interesting to observe that Eq. (|18|l remains in- 
variant if we transform S!y — > e 8 ^ SVj , where the 0j can be 
arbitrary real numbers. According to Eq. this cor- 
responds to the freedom of redefining the phases of the 
outgoing waves by different amounts on different wires. 



IV. ANALYSIS OF THE RENORMALIZATION 
GROUP EQUATIONS 

To simplify our analysis, we will make two assump- 
tions. 

(i) The couplings gu and g 2 i have the same value on all 
the wires, and therefore the subscript i on g\ and 172 can 
be dropped. 

(ii) The S'-matrix is at a FP of Eq. , so that S does 
not flow with the length scale. 

We will now consider several possibilities for the S- 
matrix, and will study the RG flows and FPs of the 
Kondo couplings Jij in each case. The different possi- 
bilities can be realized in terms of quantum wires and 
quantum dots containing the impurity spin as shown in 
Fig. 2. 




(a) 



FIG. 2: Schematic pictures of the system of wires (shown by 
solid lines), an impurity spin (shown inside a circle), and the 
coupling between the spin and the wires (dotted lines). Fig- 
ures (a) and (b) show the cases of disconnected and Griffiths 
S-matrices respectively. 



A. N disconnected wires 

The S-matrix for N disconnected wires is given by the 
A x N identity matrix (up to phases). (We will assume 
that N > 2) . A picture of the system is indicated in Fig. 
2 (a); the wires are disconnected from each other, and the 
end of each wire is coupled to the impurity spin. A more 
microscopic description of the system will be discussed 
in Sec. V. 

Let us consider a highly symmetric form of the Kondo 
coupling matrix in which all the diagonal entries are equal 
to Ji and all the off-diagonal entries are equal to J2, 
with both Ji and J2 being real. (In the language of the 



three-dimensional A-channel Kondo problem, J2 denotes 
coupling between different channels). Since the S-matrix 
is also symmetric under the exchange of any two of the 
A indices, such a symmetric form of the Kondo matrix 
will remain intact during the course of the RG flow. In 
other words, it is natural for us to choose the J matrix 
to have the same symmetry as the S'-matrix, since that 
symmetry is preserved under the RG flow. Eq. Ijl8(l gives 
the two-parameter RG equations 



dJi 
d\nL 

dJ 2 
d\nL 



1 



2itvf 
1 

2nvp 



[J 2 X +(N- 1) J 2 2 + 2g 1 J 1 



[2Ji J 2 + (A - 2) Jf 



(92 - 2 9l )J 2 ] . 

(19) 



(For A = 2 and g\ — 0, Eq. (|19(l agrees with the results 
in Ref. j^). 

Since g\ (L = 00) = 0, Eq. fli^l has only one FP at 
finite values of (Ji, J2), namely, the trivial FP at (0,0). 
We then carry out a linear stability analysis around this 
FP. [Given a RG equation of the form dX/dhiL = aX, 
we will say that the FP at X = is stable if a < 0, 
unstable if a > 0, and marginal if a = 0. In the marginal 
case, we look at the next order term; if dX/d\nL = 
bX 2 and b > 0, we say that the FP at X = is stable 
on the x < side and unstable on the x > side.] If 
v = g 2 (L = 00)/ (2itvf) > (repulsive interactions), the 
stability analysis shows that the trivial FP is stable to 
small perturbations in J2 . For small perturbations in Ji , 
this FP is marginal; a second order analysis shows that 
it is stable if Ji < and unstable if Ji > 0, i.e., it is the 
usual ferromagnetic FP which is found for Fermi liquid 
leads. However, the approach to the FP is quite different 
when the leads are TLLs. At large length scales, the FP is 
approached as Ji ~ — 1/lnL and J2 ~ 1/L V . From this, 
we can deduce the behavior at very low temperatures, 
namely, 

Jl ~ - l/ln(T K /T), and J 2 ~ {T/T K f . (20) 

where we have introduced the Kondo temperature Tk- 
(This is given as usual by Tk ~ Ae~ 27TVp / J , where A is 
an energy cut-off of the order of the Fermi energy Ep, J is 
the value of a typical Kondo coupling at the microscopic 
length scale as explained after Eq. (|22(l . and 1/(2ttUf) is 
the density of states at Ep)- The form in Eq. (|20|l is in 
contrast to the behavior of J2 for Fermi liquid leads, i.e., 
for gi = 92 = 0. In that case, Eq. I|19|l can be solved 
exactly in terms of the linear combinations Ji — J2 and 
Ji + (A - 1)J 2 ; we again find a FP at (Ji, J 2 ) = (0,0), 
with 



Ji 



l/ln(Tjc/T), and J 2 - l/\n{T K /T) 2 . (21) 



Note that J2 approaches zero faster than Ji for both 
Fermi liquid leads and TLL leads; but for the latter case, 
it goes to zero much faster, i.e., as a power of T. 
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Eq. i|21[l is valid provided that neither Ji — J 2 nor 
Ji + (N — 1) J2 is exactly equal to zero; if one of them is 
exactly zero and the other is not, then both Ji and J 2 go 
as 1/ \n(Tx/T). However, having one of the two combi- 
nations exactly equal to zero requires a special tuning in 
a microscopic model, as we will see in Sec. V. In general, 
therefore, the powers of 1/\h(Tk/T) in J\ and Ji are 
different; this does not seem to have been noted in the 
earlier literature. 




FIG. 3: RG flows of the Kondo couplings for three discon- 
nected wires, with U(0) = U(2k F ) = 0.2(2tto f ). 

Figure 3 shows a picture of the RG flows for three 
wires for U(0) = U(2k F ) = 0.2(2ttw f ). [This gives 
a value of v which is comparable to what is found in 
several experimental systems (see |45| and references 
therein). In all the pictures of RG flows, the values of 
Jij are shown in units of 2itvf-] We see that the RG 
flows take a large range of initial conditions to the FP 
at (0,0). For all other initial conditions, we see that 
there are two directions along which the Kondo couplings 
flow to large values; these are given by Jij J\ = 1 and 
J2/J1 = -l/(N - 1) (with N = 3). [On a cautionary 
note, we should remember that the RG equations stud- 
ied here are only valid at the lowest order in Jy and gi, 
i.e., for the case of weak repulsion (or attraction) and 
small Kondo couplings.] 

The fact that the Kondo couplings flow to large values 
along two particular directions can be understood as fol- 
lows. For values of J\ and Ji much larger than g\ and 
92, one can ignore the terms of order 171 and gi in Eq. 
(|19p. One then obtains the decoupled equations 



d [Ji - Ji] 1_ 

d\nL 2itvf 

d [J 1 + {N- 1) J 2 ] 1_ 

d\nL 2irvF 



(Ji - J2) 2 , 

(Ji + (N -1)J 2 ) 2 



From these equations one can deduce that the couplings 
can flow to large values in one of two ways, depending 
on the initial conditions. Either J\ + (N — l)Ji goes 
to 00 much faster than J± — Ji (this is what happens 
in the first quadrant of the figures in Figs. 3 and 4), 
or Ji — Ji goes to 00 much faster than J\ + (N — 1) Ji 
(this happens in the fourth quadrant of the figures in 
Figs. 3 and 4). A third possibility is that Ji remains 
exactly equal to zero while J\ — > 00; however, this can 
only happen if one begins with Ji exactly equal to zero. 
(This also seems to happen if the interactions are strong 
enough as we will discuss in Sec. VII). We will provide 
a physical interpretation of the first two possibilities in 
Sec. V. 

Eq. has the form dJ/dlaL = J 2 /(2nv F ). If J(d) 
denotes the value of J at a microscopic length d, and 
J{d) <C 2itvf, then it becomes of order 1 at a length 
scale Lq given by Lo/d ~ e 27™ F /j(d). corresponding 
temperature is given by Tk ~ Ae~ 27TVF ^ J ^ d \ 

Finally, note that the special case with Ji = and 
.9i = .92 = is equivalent to the Kondo problem in three 
dimensions with N channels and no coupling between 
channels j^- In the three-dimensional case, the RG 
equation has been derived to fifth order in the Kondo 
coupling 0| . This reveals a stable FP at a finite value 
of the coupling 



Ji 



N 



(23) 



where 1/(2itvf) is the density of states at the Fermi en- 
ergy. Thus the couplings Jij need not really flow to infin- 
ity as Fig. 3 would suggest; one may find strong coupling 
FPs lying at values of order 2itvf if one takes into account 
terms of higher order in the RG equations. In Sec. VII, 
we do find a strong coupling FP for sufficiently strong 
inter-electron interactions. 

Although we have discussed the case of completely dis- 
connected wires here, the results do not change signif- 
icantly if we allow a small spin-independent tunneling 
amplitude of the form 



0) ¥ jia ( Xj = 0) . (24) 



This is equivalent to changing the 5-matrix slightly away 
from the identity matrix. Using the RG equation in 114(1 , 
we find that the parameter r satisfies the RG equation 



dr 
d\nL 



1 



2itvf 



(92 - 2 5 i) T 



(25) 



(22) 



This has the same form as the interaction dependent 
terms in the RG equation for Ji in (|19|) . Hence, r also 
scales at low temperatures as T" just like J 2 in Eq. I|2U|I . 
Thus the contributions of both r and J 2 to the conduc- 
tance go as {T/T K ) 2v . 

Here and subsequently we have not discussed the case 
of attractive interactions (gi < 0) . The stability analysis 
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can easily be suitably modified in that case; some of the 
directions for the RG flows may become stable and others 
may become unstable if the sign of g 2 is reversed. 



B. Griffiths S'-matrix for N wires 



This is the case in which all the N wires are connected 
to each other and there is maximal transmission, subject 
to the constraint that there is complete symmetry be- 
tween the N wires. (We will again assume that N > 2.) 
A picture of the system is indicated in Fig. 2 (b); the 
wires are connected to each other at a junction, and the 
junction is also coupled to the impurity spin. A more 
microscopic description of the junction will be discussed 
in Sec. V. 

The maximally transmitting completely symmetric S- 
matrix is also called the Griffiths S-matrix and has all the 
diagonal entries equal to — 1+2/N and all the off-diagonal 
entries equal to 2/N. Since here, too, the S- matrix is 
fully symmetric in the N wires, we again consider the 
highly symmetric form of the Kondo coupling matrix as 
in the previous subsection, with real parameters J\ and 
J 2 as the diagonal and off-diagonal entries respectively. 
Eq. (|T%1) then gives 



dJi 
d\nL 



dJ 2 
d\nL 



JL-tf + (N-l)JS + 2. 91 d-|) 2 ^ 

-4,x(l-|) (l-^)J a ], 
^ [ 2J 1 J a + i N-2)Ji-f(l-l)J 1 



+ (g 2 -2 9l (l--) 2 ) J 2 



(26) 



(For N = 2, i.e., a full line with an impurity spin cou- 
pled to one point on the line, Eq. (|26() agrees with the 
equations derived in Ref. [2^1. 

The only FP of Eq. QXft is again the trivial FP at the 
origin. A linear stability analysis shows that this FP is 
unstable in one direction (J 2 ) and marginal in the other 
(Ji) for g 2 (L = 00) > 0. 

Figure 4 shows a picture of the RG flows for three 
wires for 17(0) = U(2k F ) = 0.2(2tto f ). We see that there 
is no stable FP at finite values of the couplings. The cou- 
plings flow to large values along one of the two directions 
J 2 /Jj = 1 and J2/J1 = -l/(N- 1). The reason for this 
is the same as that explained around Eq. I|22|l since the 
RG equations in (|19|l and 1261) have the same form for 
large values of J\ and J 2 . 




FIG. 4: RG flows of the Kondo couplings for the Griffiths 
S-matrix for three wires, with U(0) = U(2k F ) = 0.2(2nv F ). 



C. Chiral S'-matrix for three wires 



We choose a chiral S'-matrix of the form 



S 




(27) 



where 7 is a complex number satisfying (7! = 1. (We will 
see a physical realization of this form in Sec. V. Alter- 
natively, we could have considered an S'-matrix which is 
the transpose of the one given above). 

Let us consider a Kondo coupling matrix of the form 



J\ Ji J 2 
J 2 Ji J 2 
J 2 J 2 J\ 

where J\ is real but J 2 can be complex. 
Then Eq. (0 gives 



(28) 



dJx 
d\nL 

dJ 2 
d\nL 



2nvp 
1 

2nvp 



[Jl + 2\J 2 f] , 
[ 2Ji J 2 + (J 2 *) 2 



2 92J2 } 



(29) 



[Note that the above equations remain invariant under 
the transformation J 2 — > J 2 e l27T / 3 or J 2 e _l27r / 3 . We will 
see in Sec. V. C that a lattice realization of the chiral 
S-matrix has the same symmetry] 

One can again show that the only FP of Eq. I|29|l is 
the trivial FP at the origin. A linear stability analysis 
shows that the trivial FP is unstable in one direction 
(J2) and marginal in the other (Ji) for g 2 [L — 00) > 0. 
Figure 5 shows a picture of the RG flows for three wires 
for U(0) = U(2k F ) = Q.2{2izv F ). The upper and lower 
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FIG. 5: RG flows for the chiral S-matrix for three wires, with 
0(0) = U{2k F ) = 0.2(2tto f ). The upper and lower figures 
show the magnitude and phase respectively of J2. 



figures show the way in which the magnitude and phase 
of J 2 evolve. We see that there is no stable FP at finite 
values of the couplings. The phase of J 2 flows towards 
one of the three values, or ±27r/3; this is consistent 
with the symmetry of J 2 pointed out after Eq. (|29|l . 
Further, ,J\ and the magnitude of J 2 flow in such a way 
that Ji + 2 1 J 2 \ grows much faster than J\ — \Ji\. We can 
understand these observations as follows. 

For values of J\ and J2 much larger than (72, one can 
ignore the term of order g 2 in Eq. I|29(l . If we write 
J 2 = \J 2 \e l,p , we find that 



f 



dlnL 
d\J 2 \ 
d\nL 



\J 2 \ sin(3<£) , 
2Ji|J 2 | + |J 2 | 2 cos(30) ] . (30) 



The first equation in (|3(J[I shows that <f> = 0, ±7r/3, ±27r/3 
and 7r are fixed points; however, since IJ2I flows to 00 



under RG, only the values <fi = and ±27t/3 are stable. 
Substituting this fact that cos(30) — > 1 in the second 
equation in l(3()() , and combining it with the first equation 
in 1)29(1 . we obtain the decoupled equations 



d [Ji 



|J 2 



1 



d\nL 
d [Ji +2|J 2 | 



dhxL 



2itvf 
1 



i.h-\.h\? , 
(,h+2\J 2 \f 



(31) 



From this we deduce that J\ + 2\J 2 \ must flow to 00 
much faster than J\ — \J 2 \ since Ji + 2| J 2 \ > Ji — | J 2 \ 
to begin with. Note that unlike the disconnected and 
Griffiths cases, where J\ and J 2 flow to large values in 
two possible ways (with \J 2 \/Ji — ► 1 and —1/(N — 1) 
respectively), in the chiral case, J\ and J 2 flow to large 
values in only one way, along the direction \J 2 \/Ji = 1. 



INTERPRETATION IN TERMS OF LATTICE 
MODELS 



We will now see how the different iS-matrices and RG 
flows discussed in Sec. IV can be interpreted in terms 
of lattice models |2^| . This will provide us with physical 
interpretations of the various kinds of RG flows and FPs. 
We will concentrate on what the lattice models imply 
about the structure of the region near the junction, rather 
than the form of the interactions between the electrons 
in the bulk of the wires which has already been discussed 
in Sec. II. (The interactions can be introduced in the 
lattice model by, for instance, writing a Hubbard term 
at each site). We will again discuss three different cases 
here. (The models shown in Fig. 6 and discussed below 
in detail can be thought of as providing a microscopic 
picture of the systems shown in Fig. 2). 



V 



® 



3 2 10, 



(a) 



CD »---»- -< 



(b) 



FIG. 6: Lattice models for some of the S'-matrices for three 
wires, (a) can be a model for the disconnected and Griffiths 
S'-matrices, while (b) can be a model for the chiral S'-matrix. 
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A. N disconnected wires 

This system can be realized by a lattice of the form 
shown in Fig. 6 (a) . N wires meet at a junction which 
is labeled by the site number 0; all the other sites are 
labeled as n — 1,2, with n increasing as one goes 
away from the junction. (The lattice spacing will be set 
equal to one). We take the Hamiltonian to be of the 
tight-binding form, with a hopping amplitude equal to —t 
on all the bonds (where t is real), except for the bonds 
which connect the sites labeled as n = 1 on each wire 
to the junction site; we set those hopping amplitudes 
equal to zero. This is equivalent to removing the junction 
site from the system; we will therefore not consider that 
site any further in this subsection. We then obtain a 
system of disconnected wires with an S'-matrix which is 
equal to —1 times the identity matrix. To show this, 
we consider a wave which is incoming on wire i with a 
wave number k, where < k < it. We then find that 
the corresponding eigenstate of the Hamiltonian has an 
energy equal to Ek = — 2tcosfc, and a wave function 
given by 

%j) ik (n) = e~ lkn - e tkn for u=l,2, ••• on wire i , 
= at the junction and on all wires j ^= i . 

(32) 

We introduce an on-site potential which is equal to [i at 
all sites. In the absence of interactions, the ground state 
is one in which all the states with energies going from 
= — 2t up to /i are filled; the Fermi wave number hp is 
given by \x — — 2t cos kp, assuming that fi lies in the range 
[— 2t, 2t] . [We can then redefine all the wave numbers k by 
subtracting kp from them as indicated after Eq. (JTJ; the 
redefined wave numbers then run from —A to A, where 
A is of order kp.] 

Let us now consider coupling the impurity spin to the 
sites labeled as n — 1 on the different wires by the fol- 
lowing Hamiltonian 

tfspin = F 1 S- E MM) 

+ F *$- E E *£(m)ht 

(33) 

where ^ a (i, 1) denotes the second quantized electron field 
at site 1 on wire i with spin a. (Eq. (|42|l below will pro- 
vide a justification for this Hamiltonian). In Eq. <I33|) . 
i 7 ! and F 2 denote amplitudes for spin-dependent scatter- 
ing from the impurity within the same wire and between 
two different wires respectively. Namely, a spin-up elec- 
tron coming in through one wire can get scattered by the 
impurity spin as a spin-down electron either along the 
same wire (F\) or along a different wire (i 7 ^). We then 



find that the Kondo coupling matrix J^- in Eq. © is as 
follows: all the diagonal entries are given by J\ and all 
the off-diagonal entries are given by J 2 , where 

J 1 = AFi sin 2 k F , 
and J 2 = 4F 2 sin 2 fc F (34) 

for modes with redefined wave numbers lying close to 
zero. This is precisely the kind of Kondo matrix whose 
RG flows were studied in Sec. IV. A. The flows of the 
parameters J\ and J 2 considered there can be translated 
into flows of the parameters F± and F2 here. In partic- 
ular, the approach to the FP at (Ji, J2) = (0,0) given 
by Eq. (|20J) at low temperatures implies that spin-flip 
scattering within the same wire or between two different 
wires will have quite different temperature dependences. 

The flows to strong coupling shown in Fig. 3 can be 
interpreted as follows. In the first quadrant of Fig. 3, we 
see that J\ + (N — 1) J 2 goes to 00 faster than \J\ — J 2 |; 
Eq. Q34JI then implies that F\ and F 2 go to 00. In the 
fourth quadrant of Fig. 3, J\ — J 2 goes to 00 faster than 
I Ji + (N — 1) J 2 |; this implies that F\ goes to 00 and F 2 
goes to —00 as —Fi/(N — 1). 

These flows to strong coupling have the following in- 
terpretations. In the first case, F± and F2 flow to 00. 
From Eq. this implies that the impurity spin (of 

magnitude S) is strongly and antiferromagnetically cou- 
pled to only one electronic field, namely, the 'centre of 
mass' field given by ^ ■ 1)/ VN (suppressing the spin 
labels and the Pauli matrices for the moment). Hence 
that field and the impurity spin will combine to form 
an effective spin of S — 1/2. In analogy with the three- 
dimensional Kondo problem, we can say that the impu- 
rity spin is under-screened or exactly screened if S > 1/2 
or 5 = 1/2 respectively. In the second case, Fi and 
F 2 = —Fi/(N — 1) go to 00. Using Eq. (|3^|l . we can then 
show that the impurity spin is strongly and antiferromag- 
netically coupled to the N ~ 1 'difference' fields (given 
by the orthogonal combinations [^(1,1) — ^(2, l)]/y/2, 
[*(1, 1) + *(2, 1) - 2*(3, l)]/\/6, • • •)■ Hence those fields 
and the impurity spin will combine to give an effective 
spin of S-(N-l)/2 = S+l/2- JV/2. Thus the impurity 
spin is under-screened, exactly screened or over-screened 
if 25 + 1 is greater than, equal to or less than N respec- 
tively. 

B. Griffiths S'-matrix for N wires 

This system can again be realized by the lattice shown 
in Fig. 6 (a) and a tight-binding Hamiltonian. However, 
we now take the hopping amplitude to be —t on all bonds, 
except for the bonds which connect the sites labeled as 
n = 1 on each wire to the junction site; on those bonds, 
we take the hopping amplitude to be t% = — ty/2/N . 
The on-site potential is taken to be fj, at all sites, includ- 
ing the junction. We then find that the S'-matrix is of 
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the Griffiths form for all values of the wave number k. 
Namely, for a wave which is incoming on wire i with a 
wave number k, the wave function is given by 



4>ik(n) = e~ ikn - (l-^)e lfcn on wire i, 

_ ^_ & ikn Qn a jj w j reg j -J- i 

N 
2 

= — at the junction site . (35) 

We now consider coupling the impurity spin to the 
junction site labeled by zero, and the sites labeled as 
n = 1 on the different wires by the following Hamiltonian 



ffspi„ = F 3 S- £ *t (0 )£k^(o) 



a,0 



+ F *s- EE */»(«,!). 



t ct,p 



(36) 



where ^ a (0) denotes the second quantized electron field 
at the junction site with spin a. (Sec. VI will provide a 
justification for this kind of a coupling). Then the Kondo 
coupling matrix Jy in Eq. © takes the following form: 
all the diagonal entries are given by J\ and all the off- 
diagonal entries are given by J 2 , where 



4F-J 2 
Ji = -j£ + 2F 4 [1 - (l-_)cos2fc F ] 



and J 2 



4F3 

AT2 



4F 4 



cos 2kp 



(37) 



for modes with wave numbers lying close to zero. The 
RG flows of this kind of Kondo matrix were studied in 
Sec. IV. B. 

In terms of F3 and F4, the variables in Eq. (1221) are 
given by 



and Jx + {N-l)J 2 



Ji - J 2 = 2F 4 (1 - cos2fc F ) , 
N 



2F 4 (l + cos2fc F ) . 

(38) 



Since < kp < ir, l±cos2fc F lie between and 2. In the 
first quadrant of Fig. 4, J\ + ( V — 1) J 2 goes to 00 faster 
than \J\ — J 2 |; Eq. l|3*5ll then implies that F 3 goes to 00 
and \Fi\ <C F3. In the fourth quadrant of Fig. 4, Ji — J 2 
goes to 00 faster than \J\ + (N — 1) J 2 |; this implies that 
F4 goes to 00 and F3 goes to —00. 

These flows to strong coupling have the following inter- 
pretations. In the first case, F3 flows to 00 which means 
that the impurity spin (of magnitude S) is strongly and 
antiferromagnetically coupled to an electron spin at the 
junction site n = 0; hence those two spins will combine 
to form an effective spin of S — 1/2. (This case will be 
discussed in detail in Sec. VI). In the second case, F3 



goes to —00 while F4 goes to 00; hence the impurity spin 
is coupled strongly and ferromagnetically to an electron 
spin at the site n — 0, and antiferromagnetically to elec- 
tron spins at the sites labeled as n = 1 on each of the N 
wires (see Fig. 6 (a) for the site labels). Hence the im- 
purity spin will combine with those N + 1 spins to form 
an effective spin of S + 1/2 — N/2. Interestingly, we see 
that the magnitudes of the effective spins formed in the 
strong coupling limits in the first and fourth quadrants 
are the same in the cases of N disconnected wires and 
the Griffiths 5-matrix. 



C. Chiral S'-matrix for three wires 

This system can be realized by a lattice of the form 
shown in Fig. 6 (b). The three wires meet at a triangle; 
the sites on each wire are labeled as n = 1, 2, ■ • •. The 
hopping amplitude is taken to be —t on all the bonds, 
except for the three bonds on the triangle. On those 
bonds, we take the hopping amplitude to be complex, and 
of the form — te l8 in the clockwise direction and —te~ lB 
in the anticlockwise direction. [We can think of the total 
phase W of the product of hopping amplitudes around 
the triangle as being the Aharonov-Bohm phase arising 
from a magnetic flux enclosed by the triangle. Such a 
flux breaks time reversal symmetry which makes the S'- 
matrix non-symmetric. Note that since only the value of 
W modulo 2-7T has any physical significance, we are free 
to shift the value of 9 by ±27r/3. This changes the phase 
of the coupling J 2 defined below.] We then find that the 
S'-matrix is of the chiral form given by Eq. 127J1 , provided 
that the wave number k satisfies 



J(38+k) 



1 . 



(39) 



The phase 7 in (|27|l is then given by e^ 0+fc - ) . [Unlike 
the disconnected and Griffiths cases, we have not found 
a lattice model which gives an S-matrix as in (|27|) for 
all values of the wave number k.] Given a value of 9, 
we therefore choose a chemical potential /i = —2tcoskp 
such that kp satisfies Eq. (139(1 . Since the properties of 
a fermionic system at low temperatures are governed by 
the modes near kp, the above prescription produces a 
system with a chiral S-matrix. 

We now consider coupling the impurity spin to the 
three sites of the triangle through the Hamiltonian 



E E ^ mm) 

i q,/3 



(40) 

Then the Kondo coupling matrix Jij in Eq. © takes the 
form given in Eq. I|28|l . where 



(41) 



Ji = 2^5 , 
and J 2 = F 5 e-« e+3k ^ 
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for modes with wave numbers lying close to zero. This is 
a special case of the Kondo matrix given in Eq. i|28|) . [To 
obtain the most general form given in we need to 

introduce another parameter, such as a coupling of the 
impurity spin to the sites labeled by n = 2 in Fig. 6 (b) .] 
The RG flows of this kind of Kondo matrix were studied 
in Sec. IV. C. 



VI. EXPANSION AROUND A STRONG 
COUPLING FIXED POINT 

In Sec. V, we considered several examples of S- 
matrices and the RG flows of the Kondo coupling. In 
most cases, we found that the Kondo couplings flow to 
large values. We can now ask whether the vicinity of the 
strong coupling FPs can be studied in some way. We will 
see that it is possible to do so through an expansion in 
the inverse of the Kondo coupling |39| . 

We will consider one example of such an expansion 
here. Following the discussion given after Eq. let 
us assume that the RG flows for the case of the Griffiths 
S-matrix have taken us to a strong coupling FP along 
the direction J2/J1 — 1, as shown in the first quadrant 
of Fig. 4. This implies that the coupling of the impurity 
spin S to an electron spin at the junction site n = has a 
large and positive (antiferromagnetic) value Fa, while its 
coupling to the sites labeled as n = 1 on each of the wires 
has a finite value F4 which is much less than F3 (the site 
labels are shown in Fig. 6 (a)). The ground state of the 
-F3 term (namely, just the first term in Eq. (|36|0 consists 
of a single electron at site n = which forms a total spin 
of 5 — 1/2 with the impurity spin. The energy of this spin 
state is —F 3 (S + l)/2; this lies far below the high energy 
states in which there is a single electron at site n = 
which forms a total spin of 5 + 1/2 with the impurity 
spin (these states have energy F 3 S/2), or the states in 
which the site n = is empty or doubly occupied (these 
states have zero energy). 

We can now do a perturbative expansion in I/F3. We 
take the unperturbed Hamiltonian to be one in which the 
hopping amplitudes on all the bonds are —t, except for 
the bonds connecting the sites labeled as n = 1 on the 
different wires to the junction site; we take those hopping 
amplitudes to be zero. (This means that the unperturbed 
Hamiltonian corresponds to the case of N disconnected 
wires). We also include the spin coupling proportional to 
F3 in the unperturbed Hamiltonian. We take the pertur- 
bation -f/pert as consisting of (i) the hopping amplitude t\ 
on the bonds connecting the sites labeled as n = 1 to the 
junction site, and (ii) the F4 term in Eq. (J2EJ). Using this 
perturbation, we can find an effective Hamiltonian [39| . 
[Once again, we use the expression in Eq. (|15|) . where 
the high energy states are the ones listed in the previous 
paragraph. We will work up to second order in t\ and 
F4.] If 5 > 1/2, we find that the effective Hamiltonian 



has no terms of order t\ or tii 7 ^, and it is given by 

-ffeff = -Fl,cff Scff ■ ^ Si 

i 

+ c X! ( s ° s ■ ^ ef£ ■ + D / Sj 

i=£j i<j 

(42) 

plus some constants, where 

a,/3 



F hcS 
C 

and D 



t' 2 



(5 + 1) (25+1) F 3 
2 (S + 1) Ft 



2(S + 1)F 4 
25+1 



(25 + 1)3 F 3 



(5 + 1) (25+1) ^3 
2Fi 



(25 + l) 3 F 3 ' 

n 

(25 + 1) F 3 



(43) 



In Eq. H42JI. S c s denotes an object with spin 5 — 1/2. 
We thus find a weak interaction between the spin 5—1/2 
and all the sites which are nearest neighbors of the site 
n = as shown in Fig. 6 (a). 

If the impurity is a spin- 1/2 object (i.e., 5 = 1/2), 
then the electron at the site n = forms a singlet with 
the impurity. In that case, only the last term in Eq. 
(|42H survives. However, there are other terms in the ef- 
fective Hamiltonian which are of higher order in ti /F3 
than in (|42() ; these have been calculated in Ref. |47| for 
the case 5 = 1/2. One of these terms describes spin- 
independent tunneling from one wire to another, of the 

form Ei/jEa*!!*' 1 )*"^ 1 )' Tnis is a contribution 
to the 5-matrix at the junction, and it can contribute 
to the conductance from one wire to another as we will 
discuss in Sec. VIII. 

Returning to the case 5 > 1/2, we note that the last 
two terms in Eq. (|42|) are irrelevant as boundary opera- 
tors if 52 {L — oo)/(2itvf) is small; this is because Si has 
the scaling dimension 1 — §2/ {2'kvf) (as one can see from 
Eq. H19I) - ). and therefore the product Sj ® Sj has the scal- 
ing dimension 2(1 — {litvp)) which is larger than 1. 
The first two terms in Eq. I|42|) have the same form as in 
Eqs. and (|34|l . where the effective Kondo couplings 



Ji, e f£ = 4Fi :Cff sin 2 k F 
and J2,eS = 4F 2 ^ c ff sin 2 fcj? 



(44) 



are equal, negative and small. We can now study the RG 
flow of this as in Sec. IV. A. With these initial conditions, 
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Eq. (|19fl and Fig. 3 show that the Kondo couplings flow 
to the FP at (Ji ieff , J 2 .es) = (0, 0). 

In this example, therefore, we obtain a picture of the 
RG flows at both short and large length scales. We start 
with the Griffiths S-matrix with certain values of the 
Kondo coupling matrix, and we eventually end at the 
stable FP of the disconnected S-matrix for repulsive in- 
teractions, g 2 {L — oo) > 0. 

We will not discuss here what happens for the other 
possible RG flow for the Griffiths S'-matrix, in which 
Ji and J 2 become large along the direction J 2 /J\ = 
-l/(N - 1). As we noted in Sec. V B, N + 1 spins get 
coupled strongly to the impurity spin in that case; an 
expansion in the inverse coupling is much more involved 
in that case. For the same reason, we will not discuss 
expansions in the inverse coupling for the flows to strong 
coupling for the disconnected and chiral S-matrices. 



VII. DECOUPLED WIRES WITH STRONG 
INTERACTIONS 

In this section, we will briefly discuss what happens 
if there are N decoupled wires and the interactions are 
strong. For the decoupled S'-matrix, one can 'unfold' the 
electron field in each semi-infinite wire to obtain a chi- 
ral electron field in an infinite wire, and then bosonizc 
that chiral field 0, 0, H3] • In the language of bosoniza- 
tion, the interaction parameters are given by K p for the 
charge sector and K a for the spin sector. Spin rotation 
invariance implies that K a = 1, while K p is related to 
our parameters gi as follows |2(ij |. 



K n = 




g\ I Trvp + (gi - 2g 2 )/2irv F 



1 



g^j-KVF 

gi - 2g 2 
2ttvf 



(gi - 2g 2 )/2-Kv F 



(45) 



In the second line of the above equation, we have taken 
the limit of small gi since we have worked to lowest order 
in the gi in the earlier sections. From Eq. (|llf> . we see 
that 2g 2 — <?i is invariant under the RG flow. The case 
of repulsive interactions corresponds to 2g 2 ~ g\ > 0, i.e., 
K p < 1. 

The case of two decoupled wires (N = 2) has been 
studied by Fabrizio and Gogolin in Ref. |2flj. They 
showed that if the interactions are weak enough, the 
Kondo couplings J± and J 2 are both relevant; their re- 
sults then agree with those discussed in Sec. IV A. But 
if the interactions are sufficiently strong, i.e., K p < 1/2, 
then 3 2 is irrelevant and flows to zero. 

We will now show that their results can be generalized 
to the case of N wires; one finds that there is again a 
value of K p below which J 2 is irrelevant. Following Ref. 
[stij . we can write the spin-up and down Fermi fields ^i : u 
in wire i in terms of the charge and spin bosonic fields 



$i iP and ^i :a . Close to the junction denoted as Xj 



we have 



and * 



, (46) 



where we have used the fact that K a = 1, and we have 
not explicitly written the arguments of the fields {xi = 0) 
for notational convenience. The r\i^ a denote Klein factors, 
and d is a short distance cut-off; these will not play any 
role below. 

In bosonic language, the Hamiltonian H — H + H lnt 
in Eqs. (@J and 10 is given by 



H = 



-E 

4tt ^ 



dxi 



dxi 



V dxi 



. (47) 

where v p , v a denote the charge and spin velocities re- 
spectively. The bosonic fields satisfy the commutation 
relations 



d$i,g{Xi) 

dx,. 



, ®j b(xj) ] = i 2-7T S (l b Sij S(xi - Xj) , (48) 



where a,b = p. a. 

The impurity spin part of the Hamiltonian is given by 



E E *L%^ 



i,0 



E E 



3,P 



(49) 



The spin densities on different wires are given by 

- f tft + tf-t - * f i 1 - 1 dq>l ' a 
2 [ J 2V27T dxi ■ 

The other terms take the form 

U/t xb- t ^ P (»/V2)[-*., P /\/^- <i> ^+' !> i,p/v / ^+' ! '^] 



(50) 



,(i/V2)[-*i,p/\/K7-*i,<r+*i.p/V^p-*^«'] 



(51) 



and so on. In 149(1 and 1(51(1 . we have not explictly writ- 
ten the arguments of the fields, Xi = Xj = 0; we will 
continue to do this wherever convenient. [The bosonic 
forms of the fermion bilinears in Eqs. ((501) and ((51(1 are 
so different because we are using abelian bosonization. 
For the same reason, we will find it useful to distinguish 
between the different components of J\ and J 2 , i.e., Ji 2 , 
Jix, Jiz, and J 2 ±.] Let us define N 'orthonormal' linear 
combinations of the spin boson fields, namely, the 'centre 
of mass' combination 



(52) 
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and the 'difference' fields 



n 

* s = , 1 , n [ E *w 



n $„ 



+l,a 



(53) 



where n = 1, 2, • • • , TV — 1. We can now write Eq. I|49l) in 
the bosonic language. We obtain 



H, 



Nm ' " 2V27T 

Jl 

And 



Ju sz y d * 



dx. 



Ju_ [s+ e ,vw*° y e iEn«? 

4W L 



7rd 



5* £ sin (E 6 « $ ") 

i<j \ n / 



x sin 



^,P - *j,P 



2nd 



E 



e ^» «3 ■ " cos ' 



(54) 



where the sums over n in the second, third and last lines 
run over the 'difference' fields The constants a", b 7 ^ 
and c™ in Eq. 1)54(1 satisfy the relations 

E «) 2 = 2 - | ' 

n 

E ( & S) 2 = i . 



and y {c-f = 1 - 



TV 



(55) 



for all values o(i,j. 

We can remove the phase factors exp(iy2/TV<I>°) in 
Eq. I|54[) by performing an unitary transformation of the 
total Hamiltonian H to t given by the sum of 1(47(1 and (|49ll , 
namely, H tot UH tot W 0, where 



cos 



H.c. 



(57) 



where A = Ji z — Anv a /N. We can now study the problem 
in the vicinity of the point A = Ji± = J22 = J2A. = 0. 
Note that this is a strong coupling FP, since A = implies 
that 



Ju 



4nV a 

N 



(58) 



Since the scaling dimension of e l/3 * i a is given by /3 2 /2, 
for a = p, a, we see from Eq. (I55|l that the operators 
multiplying Jix, J2z and fo± in Eq. 157(1 have the scaling 
dimensions 1-l/N, 1/2 + l/(2K p ) and 1/2 - 1/TV + 
1/{2K P ) respectively. This impies that the J\± operator 
is always relevant, while the Ji z operator is irrelevant 
if K p < 1 (repulsive interactions). Most interestingly, 
the J21. operator is relevant or irrelevant depending on 
whether K p > or < N/(N + 2). For TV = 2, this gives the 
critical value of K p to be 1/2 j^O], while for TV — * 00, the 
critical value of K p approaches 1, i.e., the limit of weak 
repulsive interactions. 

We saw in Sec. IV A that a flow to strong coupling 
is indeed possible along the line J2 = 0, although that 
line is unstable to small perturbations in J2 . We now see 
that the line is stabilized (to first order in the couplings) 
if the interactions are sufficiently strong, i.e., if 



TV 



(59) 



If J 2 flows to zero and J\ flows to large values, Eq. j3^|) 
shows that the impurity spin is coupled strongly and anti- 
ferromagnetically to the electron fields ^f(i,l) on all the 
TV wires; hence they will combine to form an effective 
spin of S — N/2. (If S < N/2, the impurity spin is over- 
screened). This describes a TV-channel AFM FP with no 
coupling between channels jU [2^ . 



VIII. CONDUCTANCE CALCULATIONS 



After this transformation, Eq. ((54(1 takes the form 



H. 



And 

J27. 
nd 



A gz yd* 



dxi 



[* + E 



H.c. ] 



E sin fe^) 

i<j \ n / 



X sm 



Our calculations for the Kondo couplings can be ex- 
plicitly applied to various geometries of quantum wires 
and a quantum dot (containing the impurity spin) shown 
in Fig. 2, such as (a) a dot coupled independently to 
each wire (disconnected S-matrix for the wires), so that 
the conductance can only occur through the dot, or (b) a 
side-coupled dot (Griffiths S'-matrix for the wires) , where 
the conductance can occur directly between the wires. In 
general, of course, one can have any S'-matrix at the junc- 
tion, so that the conductance can occur both through the 
dot and directly between the wires. 

Let us now consider the conductance near the different 
FPs [24I l35j | for the case of weak interactions. In the 
Griffiths case where the conductance can occur directly 
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between the wires, let us consider the case of small values 
of Ji, J 2 (both much smaller than 2ttvf), and g 2 3> g\. 
At high temperatures, before the J^'s have grown very 
much under RG, we see from Eq. (|26[) that J\ remains 
small, while J 2 grows due to the term g 2 J 2 . Namely, 
J 2 - (T/T K )- U , where v = g 2 (L = oc)/(2irv F ). The 
effect of J 2 is to scatter the electrons from the impurity 
spin, and thereby reduce the conductance between any 
two wires from the maximal value of Go = (4/iV 2 )e 2 /h. 
Since the scattering probability is proportional to j| , the 
conductance at high temperatures (T 3> Tk) is given by 



G - Go Go 5(5 + 1) (T/T K y 



(60) 



[The factor of 5(5 + 1) appears for the following reason. 
Consider an electron coming in through wire «; it can 
have spin up or down, and the impurity spin can have any 
value of S z from 5 to —5. We assign all these 2(25 + 1) 
states the same probability. As a result of the Kondo 
coupling J2, the electron can scatter to a different wire 
j; as a result, its spin may or may not flip, and the value 
of S z for the impurity spin can also change by or ±1. If 
we calculate the probabilities of all the different possible 
processes and add them, we get a factor of 5(5 + 1)]. 
Using Eq. (|4*5)l . we see that (fHOf) takes the form 



G — Gq 



Go 5(5 + 1) (T/T K ) 



K„-l 



(61) 



where we have used the RG flow to set g\ = and v = 
g 2 (L = 00)/ (2ttvf). On the other hand, if the leads were 
Fermi liquids (g\ = g 2 — 0), J 2 would be given by Eq. 
(I2T)| . and we would get 



G — Gq 



G 5(5 + 1) 
\n{T/T K y 



(62) 



At low temperatures, the Kondo couplings flow to large 
values; as discussed at the end of Sec. VI, their behaviors 
are then governed by the FP at ( Ji iC ff, J 2 , off) = (0, 0) of 
the disconnected wire case with an effective spin 5 G ff — 
5 — 1/2. In this case, only Jf contributes to the conduc- 
tance between two different wires. From Eq. (|2U|I . we see 
that the conductance is given by 



G ~ Go 5 off (5 off + 1) {T/T K ) 2v 
~ Go 5 off (5 off + l) {T/T K ) 1 - K ^ 



(63) 



for T <C Tk- For Fermi liquid leads, Eq. (|2~T)l implies 
that the conductance is given by 



G 



Go5 c ff (5 c ff + 1) 
HT/T K ) 4 



(64) 



Thus a measurement of the temperature dependence of 
the conductance should be able to distinguish between 
the Fermi liquid and TLL cases at both high and low 
temperatures. For the case N — 2, the expressions in 
Eqs. 1)61(1 and H63[l agree with those given in Refs. [2^, 



|35| . but Eqs. (|62|l and l|64|) differ from the expressions 
given in earlier papers (like Ref. (35|) for the powers of 
l/ln(T/T K )- (As we had discussed earlier after Eq. J2U, 
we would get the same powers of l/ln(T/Tx) as in Ref. 
[35} if J 2 was exactly equal to J\ or —Ji/(N — 1)). 

The above expressions for the conductance shows that 
for both Fermi liquid leads and TLL leads (with repul- 
sive interactions), and for both T ^> Tk and T <C Tk, 
the conductance increases with the temperature. It is 
then natural to assume that this would be true for inter- 
mediate temperatures as well, so that the conductance 
increases monotonically with temperature from to Go; 
this would be consistent with the results in Refs. [iflLlsij . 
It may be useful to discuss here why there is no Kondo 
resonance peak in the conductance at low temperatures 
in our model, in contrast to what is found in other mod- 
els (for instance Refs. 0,ESH3) an d observed experi- 
mentally [fill 152]. In our model, once the impurity spin 
gets very strongly coupled to the junction site in Fig. 
6 (b) (due to the flow to large J\ and J 2 in the Grif- 
fiths case), that site decouples from the wires; this leaves 
no other pathway for the electrons to transmit from one 
wire to another. In contrast to this, if the junction region 
was more complicated (for instance, if there were addi- 
tional bonds which connect different wires without going 
through the impurity spin, or there was a dot with several 
energy levels through which the electron can transmit), 
then the electron may still be able to transmit even after 
the impurity quenches the electron on a single site (or 
level). Hence, it may be possible for the conductance to 
increase to the unitarity limit at the lowest temperatures; 
this is known to occur for models with Fermi liquid leads. 
For TLL leads, however, our analysis remains valid even 
if there are additional bonds between the wires, because 
any such direct tunneling amplitudes are irrelevant and 
renormalize to zero as shown in Eq. Q25JI . 

Finally, let us briefly consider the case of strong inter- 
electron interactions. For K p < N/(N + 2), we saw in 
Sec. VII that a multi-channel FP gets stabilized in the 
case of N disconnected wires. To obtain the conductance 
at this point, we need to study the operators perturbing 
this point, similar to the analysis in Refs. j2f| 0, 
this has not yet been done. 



IX. CONCLUSIONS 

To summarize our results, we have studied systems of 
TLL wires which meet at a junction. The junction is 
described by a spin-independent S'-matrix, and there is 
an impurity spin which is coupled isotropically to the 
electrons in the neighborhood of the junction. The 5- 
matrix and the Kondo coupling matrix Jy satisfy certain 
RG equations. We have studied the RG flows of the 
Kondo couplings for a variety of FPs of the iS-matrices. 
Although the Kondo couplings generally grow large, one 
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can sometimes study the system through an expansion in 
the inverse of the coupling. This leads to a new system 
in which the effective Kondo couplings are weak; the RG 
flows of these effective couplings can then be studied. 

For example, at the fully connected or Griffiths S- 
matrix, we find that for a range of initial conditions, the 
Kondo couplings can flow to a strong coupling FP along 
the direction J2/J1 = 1, where their fate is decided by 
a 1/J analysis. This analysis then shows that the cou- 
plings flow to the FM FP of the disconnected S-matrix 
lying at ( Ji, e ff , Ja.eff) = (0, 0). For this system, therefore, 
one obtains a description of the system at both short and 
large length scales. 

For the case of disconnected wires and repulsive in- 
teractions, there is a range of Kondo couplings which 
flow towards a multi-channel FM FP at ( Ji, J 2 ) = (0, 0). 
At low temperatures, we find spin-flip scattering pro- 
cesses with temperature dependences which are dictated 
by both the Kondo effect and the inter-electron interac- 
tions. It may be possible to observe such scatterings by 
placing a quantum dot with a spin at a junction of several 
wires with interacting electrons. 

For other initial conditions for the disconnected case, 
the Kondo couplings flow towards the strong coupling 
FPs at Ji, I J2I — > 00. In general, this is just the single 
channel strong coupling AFM FP. But there is a spe- 
cial line where J\ — > 00 and J2 = 0; this is the multi- 
channel AFM FP. The RG equations show that both Ji 
and J2 are relevant around the weak coupling FP if the 
interactions are weak. However, if the interactions are 
sufficiently strong (i.e., K p < N/(N + 2)), we find that 
J2 — *■ 0, and the multi-channel FP gets stabilized. 

Experiments are underway to look for multi-channel 
FPs, and proposals have been made for minimizing the 
couplings between channels using gate voltages [2^]. We 
suggest here that enhancing inter-electron interactions in 
the wires offers another way of reducing the inter-channel 
coupling and thereby observing the effects of the multi- 
channel FP. 

Finally, we have discussed the temperature depen- 
dences of the conductances close to the disconnected and 
Griffiths S'-matrices, and showed that this also provides 
a way to distinguish between Fermi liquid leads and TLL 
leads. 
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